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The large iV limit of fermionic vectors models is studied using bilocal variables, in the 
framework of a collective field theory approach. The large iV configuration is determined 
completely using only classical solutions of the model. Further, the Bethe-Salpeter equa- 
tions of the model are cast as a Green's function problem. One of the main results of this 
work is to show that this Green's function is in fact the large iV bilocal itself. 
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1. Introduction 

Systematic large iV expansions remain among the most promising tools that can be 
used to probe the non-perturbative aspects of quantum field theories. For matrix theories, 
the leading contribution is determined by summing all of the planar diagrams. Except for 
the simplest models, this is still not possible. In contrast, the situation for vector models is 
much better (some work on vector /matrix models includes, for instance, [1]]). Indeed, the 
large N limit can often be constructed explicitely. The AdS/CFT correspondence relates 
the large iV limit of matrix theories to a dual gravitational theory. One hopes, that with 
a better understanding of this correspondence, large iV expansions of matrix theories will 
become possible. It would be useful to have a simpler example with which the correspon- 
dence can be studied. Recently, it has been suggested that the large iV expansions of 
vector models can also be related to dual gravitational theories @, the so-called higher 
spin theories. In much the same way that vector models were used to understand many of 
the features expected in QCD, it may be possible to use the vector model/higher spin the- 
ories correspondence as toy models for the AdS / CFT correspondence. Given this renewed 
interest in the large iV limit of vector models, it is appropriate to revisit the outstanding 
issues in the large N expansions of these theories. 

A concrete and very interesting proposal for a derivation of the vector model/higher 
spin theories correspondence was given in || . This article employed collective field theory 
to give a description of bosonic vector models in terms of bilocal variables. The collective 
field theory is a promising tool to attack this problem. Experience with the AdS/CFT 
correspondence suggests that supersymmetric examples may be the most amenable to 
analysis. For this reason, the construction of the collective field theory for fermionic vector 
models is well motivated. In particular, the descriptions of solitonic objects in the large N 
limit of fermionic vector models is not well understood. This was the primary motivation 
for this work. 

We discuss the description of the soliton of the Gross-Neveu model, both at the lead- 
ing and next to leading orders. Suprisingly, both are determined largely by the classical 
solutions of the model. This is likely to provide useful insights into how one would ap- 
proach similar questions for supersymmetric vector models which are relevant for the vector 
model/higher spin theories correspondence. Quite apart from the motivation provided by 
this correspondence, the fluctuations at sub leading order are interesting in their own right. 
Indeed, they supply important information about the stability of the solution (absence or 
prescence of tachyons) and about the space of solutions (via the zero modes) . 



1 



This paper is organized as follows: in section 2 the bilocal formalism is discussed 
and some central results stated. Using the bilocal formalism, a leading order Ansatz is 
constructed in section 3, first for the perturbative vacuum of the Gross-Neveu model, and 
then for the soliton. Section 4 is concerned with the next to leading order fluctuations. In 
section 5, we discuss our results. 

2. The bilocal formalism 

We follow an approach based on collective field theory |4[] to treat the large N limit 
of vector models ||, ||. In this approach a change in variables is implemented from the 
original fermionic field to time-ordered bilocals. The bilocal fields are singlets so that 
in terms of these variables, all iV dependence is explicit. Consequently, these are the 
correct variables in which to develop a systematic 1/N expansion. As usual, this change 
of variables gives rise to a Jacobian [|j . Incorporating the effects of this Jacobian produces 
an effective action that can be used to generate systematic corrections in 1/N. 

We start with an action invariant under U(N) transformations. These transformations 
act as (a is a spinor index) 

i> a ^il>' a = Ui> a , ^ a ^^ = ^ a U\ UeU(N) 

on our fermionic variables. We change variables from the original fermionc fields to the 
invariant bilocals, given by 

Requiring that the Schwinger-Dyson equations for arbitrary singlet operators derived using 
the bilocal variables agree with the equations derived using the orginial variables yields a 
differential equation for the Jacobian J. The equation can be solved exactly to obtain || 

log J = - [N + mL d 5 d (0)] Tr log a, 

where the trace runs both over the spinor indices and over functional space. In the above 
expression, m is the dimension of the Clifford algebra, L d is the volume of spacetime and 
the delta function is in momentum space. The leading term in this expression was also 
obtained in |7]]. In terms of this Jacobian the effective action is given by 
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S eff = -i\ogJ + NS 

= ziVTr log a + NS + imL d 5 d (0)Trlog a 
= NS + S 1 . 

It was shown in [BJ that it is consistent to treat the last term as subleading in a systematic 
1/N expansion. In the above, the fields have been rescaled such that an overall factor N 
multiplies the action. Consequently, the leading large iV configuration, a , is now obtained 
by solving 



SSq 



da 



0. 



The systematic expansion is now obtained by expanding about this leading configuration 



<r a p{x, V) = V) + ~7=Vap(x, y). 



Performing a systematic expansion of the effective action we obtain 



S eff = NS (a°) + S 1 (a°) - % -C 2 + \d 2 



oo 1 



n=l 



i(-l) 



n+1 



; ^- + mL a d(0) — 

(n + 2) y ! n 



+ 



D 



n+2 



(n + 2)! 



f2.1) 



where 



D r 



d d Xl ■ ■■ d d x n / d d yi ■■■ d d y 



5 n S 



' Sa^fc (xi, y x ) ■ ■ ■ Sa an p n (x n , y n ) 



and 



C n = Tr((a°)- 1 V ) n . 

The first corrections to the leading configuration will come from the (quadratic in 77) terms 
C<i and D2. The use of these terms to generate corrections is one of the questions with 
which we will be concerned in this paper. In particular, we solve the Bethe-Salpeter which 
is obtained as the wave equation for the bilocal following from the quadratic action. 
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3. Large N and the classical solution 

The results of the previous section (in particular, equation (|2.1|) ) clearly show that the 
effective action for the bilocal field is determined by the leading large iV configuration a . 
In this section, our aim is to argue that this leading configuration can be constructed using 
the classical solutions of the fermionic vector model. This extends earlier results obtained 
for bosonic vector models @, ||. The result for fermionic models is based largely on the 
unpublished works [fl, 110]. It extends and completes the discussion appearing in [ITT . 



As we mentioned in the introduction, our arguments are illsutrated using the Gross- 
Neveu (GN) model. However, we would like to emphasize that the methods developed in 
this article are general. The Lagrangian of the GN model is given by |12[ 



l = r( x ) {vfdj r(x) + l - g 2 (r( x )r(x)) 2 , 

where a = 1, N is a color index. Note the presence of a quartic interaction and absence 
of a mass term. In the context of a path-integral quantization, one would treat if) and ip as 
Grassman valued fields. It is thus natural to suspect that the relevant classical solutions 
are obtained by treating ijj as a classical Grassman field. We will argue that this is not the 
case - one should treat ijj as a normal commuting function. 

Recall that for a vector model, the large iV limit is reached by holding A = g 2 N fixed 
and taking N — ► oo. The leading bilocal configuration is determined by the Schwinger- 
Dyson equations. The simplification of the theory in the large N limit is reflected in the 
factorization of expectation values of singlets. The Schwinger-Dyson equations 



for the GN model are 



N5 aP 5 2 (x - y) (1) = i (^(!/)i(7T 7 ^;W+^(!/)C(^;(^!(i)) • 
In terms of the bilocal 

a a p(x,y) =C(x)^(y), 
the large iV Schwinger-Dyson equationB reads 

1 That is, the Schwinger-Dyson equation obtained after using factorization in the large N limit. 
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= iN5 a p5 2 (x -y)+i (^) au d xp op u (y, x) + g 2 o fioL {y, x)a pp (x, x), 

where the partial derivative is acting on x. To interpret this equation, note that the 
classical equation of motion reads 

= i (YT" d XfJ ,a/3v(y, x) + g 2 <Jp a (y, x)a pp (x, x). 

Thus, the large N Schwinger-Dyson equation gives the classical dynamics in the presence 
of a point source. We would like to emphasize that this result is obtained by treating ip as 
an ordinary commuting function and not as a Grassman valued field. Further, the point 
source does not arise from the specific action used. It is a general result that the large N 
configuration of the bilocal singlet field is simply the classical bilocal singlet field, obtained 
in the presence of a point source with strength proportional to N 0. 

We now turn to the problem of explicitely determining the leading bilocal configuration 
in terms of the classical solutions. It is simplest to illustrate our method on the perturbative 
vacuum (no soliton) of the theory. Towards this end, note that there are two linearly 
independent classical solutions for any given energy 

where 

u i = k i + (fi ,2 c r aa(a^, x)) 2 = k 2 + m 2 . 
Using these two solutions, it is natural to make the Ansatz 

*ap(x, y) = 9{t x - t y ) + m) ap ^-^(t,-*,)-^-*)) 
- 6{t v - t x ) ( 7 % - m) ap 
It is clear that, for unequal times, this Ansatz solves the SD equation, since then the SD 
equation reduces to the classical equation of motion. 

Integrating the SD equation with respect to t x form t y — e to t y + e 



/ dt x N8 a p8{x-y)8(t x -t % 

Jt,, — s 



t y -e 

V+S dt x [(i(-r°) av d ta - iiYT'd^apviy, x) + g 2 o Pa {y, x)o pp {x, x)] . 



t v -e 



The expression on the left is an integral over a delta function and is trivial to perform. 
Consider now the right side. It is only the first term in the above expression, the time 



derivative, which contributes. We obtain 



-iN8 af3 8(x -y) = itf) av <rp v {y t x) |££ 

i 

-0(t y - t x )i{l T v (7% - ™)p v Ale^v-'V 
This last expression is rearranged to become 



-iN6 a p6(x-y) = 25 a p^A 2 i u i e 1 



(ki(y-x)) 



which can be solved to give 



At 



N 
2Lujj 



Substitution of uf = kf + (<? 2 cr QQ (x, x)) 2 and the normalization condition ^7 
the constraint 



1 gives 
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L ° i V k i + {g 2 o aa {x,x)) 2 

which is the gap equation for this theory |J . Although we have considered the perturbative 
vacuum of the field theory the argument is general. Indeed, a completely parallel argument 
allows one to construct the large N bilocal corresponding to the kink. Using the classical 



solution of |13|] we easily obtain 



vpoiy, x) = ^2[N9(t x - t y )ipp(y, k)ip e {x, k) - N6(t y - t x )ip cf3 (y, k)ip ce (x, fc)] s 



where ijjg(x, k) is the classical solution of 



V 1 (x)+iV 2 (x) 
U l {x)+iU 2 {x) 
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tanh(mx) + 



(1 + tanh(mx)) 
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1 - 



(1 + tanh(mx)) 



3 i(cJfet— kx) 
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and ip c is related to ifi by charge conjugation. One can once again verify that this bilocal 
solves the large N SD equation ||, []TIJ. Notice that this is in fact the solution to the 



classical Gross-Neveu model equations with N = 1. However, in the construction of the 
bilocal, it is clear that the sum over the spectral parameter k is playing the role of a 
sum over color. This is reminiscent of the quenched momentum prescription for matrix 
models[[n]. Exactly the same behaviour was noticed even earlier in the context of bosonic 
vector models [fj, §|. 



4. Fluctuations about the Leading Kink Solution 

In the previous two sections we have provided a construction for the bilocal repre- 
senting the leading large iV configuration. At finite N, corrections to the leading bilocal 
will become important. In this section, we address the problem of computing these cor- 
rections. The leading correction is determined by the terms in the effective action that 
are quadratic in r\. These terms, when varied with respect to r\ give the Bethe-Salpeter 
equation of the theory, in the one soliton background. The result of this section is that 
we are able to solve this Bethe-Salpeter equation. Once again, all that is needed for this 
construction is a knowledge of the classical solutions of the theory! It is primarily with 
the analysis of this section in mind that we chose to study the Gross-Neveu model. Using 
the fact that the model is integrable, we will be able to argue that we have constrcuted 
a complete set of solutions to the Bethe-Salpeter equation. Although our construction is 
completely general!!, the proof of the completeness of the solution set makes crucial use of 
the integrability of the model. 

Our first observation is that the Bethe-Salpeter equation follows by minimizing the 
quadratic action or by linearizing the Schwinger-Dyson equations. Our use of the integra- 
bility of the Gross-Neveu model is to use the scattering data as suitable parameters for 
this linearization. Since the scattering data are a complete set, we are sure to obtain a 
complete set of solutions to the Bethe-Salpeter equation. Our approach is motivated by 
the stability analysis of |0J. 

The scalar field 

2 That is, as long as the leading configuration is known, our construction can be carried 
out. Integrability of the model is not used; once the leading large iV configuration is known the 
construction is (in principle) mechanical. 
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cr(x, t) = ipip, 
can be used to define a new field u(x, t) as 



a(x, i) 



i(x,t) 



Using the equations of motion for the Gross- Neveu model, one finds that u(x, t) satisfies 
the Sine-Gordon equation. It is possible to explicitely obtain the dependence of u(x, t) 
(and hence a(x, £)) on the scattering data, using the methods developed in []I6| . The field 
cr(x) appears in the large N Schwinger-Dyson equations as a "potential" 



= iN5 a p5 2 (x -y) + i (tT" d XfI ap u (y, x) + g 2 ap a (y, x)a(x). 

Denote the generic scattering data parameter by A. Under A — > A + AA, we obtain a 
second solution to the Schwinger-Dyson equations. This follows as a consequence of the 
fact that the Schwinger-Dyson equations have no expicit dependence on A. It is a simple 
task to show, by expanding, 

a a p(x, y;A + AA) = a a p(x, y; A) H ^ AA, 

and using the fact that both cr Q( g(x, y; A+AA) and a a p(x, y; A) satisfy the Schwinger-Dyson 
equations, that 5lTa i' 3 ^ v '^ satisfies the Bethe-Salpeter equation. Now, for the problem 
considered here, since we know the dependence of cr(x) on A (thanks to the connection 
to the Sine-Gordon model) we can compute the variation of ip^(x) using the (linearized) 
equation of motion 



(d t + d x )5i(j2 = -g 2 (a5ifj 1 + 8aip x ) 
(d t -d x )5i; 1 =g 2 (a5if; 2 + 5ai; 2 ). 
The variation of the bilocal then follows 



(4.1) 



Sa aP (x,y) = 5r a (x)r P (y)+r a (x)5r P (y). 



The solution to the system of equations (|4.1|) can be written as 



5ipi(x) 
5ijj 2 (x) 



d 2 y 



Gn(x,y) G 12 (x,y) 
G 2 \{x,y) G 22 {x,y) 



8^ 2 {y) 
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where the Green's function G oc /3(x,y) satisfies 

H^Gix, y) + g 2 a(x)G{x, y) = 5^(x - y). 

The central observation of this section is that the above equation is nothing but the 
equation for the large N bilocal and as a result, the above Green's function is equal 
to the large iV bilocal configuration. 

For the case of the Gross-Neveu model, it is possible to explicitely implement this 
construction. Further, the construction of a complete set of solutions to the Bethe-Salpeter 
equation can be used to compute the propagator for the fluctuation r\ in the background 
of the soliton. This allows one to compute (for example) the leading 1/N correction to 



the soliton mass which is known H17| . The specific form of the solutions to the Bethe- 



Salpeter equation, the use of these equation to determine the 77 propagator and corrections 



computed using this propagator will be given in fl8|j . 

If the model was not integrable, the form of 5a would have been unknown. In this 
case the large iV bilocal configuration continues to furnish a suitable Green's function. 
Using this Green's function one obtains an integral equation that can be iterated in much 
the same way that one treats scattering problems in non-relativistic quantum mechanics. 
In this case, one obtains the analog of a Born approximation for the solutions of the 
Bethe-Salpeter equation. 



5. Conclusion 

We have studied the large N limit of fermionic vectors models using a collective field 
theory approach. The large iV configuration was determined completely using only classical 
solutions of the equations of motion, obtained by treating the fermionic fields as ordinary 
commuting functions. For the case of the Gross-Neveu model, we have demonstrated that 
a complete set of solutions to the Bethe-Salpeter equation can be constructed. The central 
insight is that the solutions can be determined using a Green's function method, and 
further that this Green's function is in fact the large iV bilocal configuration itself! For 
more general theories, which will not be integrable, this insight will provide an iterative 
approximation scheme for the solutions to the Bethe-Salpeter equation. It is surprising 
that so much of the dynamics of a fermionic quantum field theory can be captured by well 
chosen classical solutions. 
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Our analysis of the Gross-Neveu model may now be extended to construct the propa- 
gator for the soliton, and compute, for example, corrections to quantities of interest which 



can be read from the (known) exact S'-matrix of the model [I7|. This work is in progress 
||18|| . It would be interesting to extend this analysis to other completely integrable models. 

For models which are not integrable, we have outlined a procedure to determine an 
integral equation for the solutions to the Bethe-Salpeter quation. The kernel of this in- 
tegral equation is the large iV bilocal. This integral equation could be iterated to obtain 
approximate solutions. The accuracy and convergence of this iterative procedure is an 
interesting question which deserves to be explored. 

I would like to express my gratitude to R. de Mello Koch and J. A. P. Rodrigues for 
discussion aiding in the completion of this work. 
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